Transmission Character of General Function Photonic Crystals 
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In the paper, we present a new general function photonic crystals (GFPCs), which refractive index 
of medium is a arbitrary function of space position. Unlike conventional photonic crystals (PCs), 
which structure grow from two mediums A and B, with different constant refractive indexes n a and 
nij. Based on Fermat principle, we give the motion equations of light in one-dimensional GFPCs, and 
calculate its transfer matrix, which is different from the conventional PCs. We choose the linearity 
refractive index function for two mediums A and B, and find the transmissivity of one-dimensional 
GFPCs can be much larger or smaller than 1 for different slope linearity refractive index function, 
which is different from the transmissivity of conventional PCs (its transmissivity is in the range of 
and 1). Otherwise, we study the effect of different incident angles, the number of periods and opti- 
cal thickness on the transmissivity, and obtain some new results different from the conventional PCs. 
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1. Introduction 

Photonic crystals (PCs) are composite structures with 
a periodic arrangement of materials with different refrac- 
tive indices in one-dimension (ID), two-dimension (2D) 
or three-dimension (3D). Due to the introduced periodic- 
ity, multiple Bragg scatterings from each unit cell may 
open a photonic band gaps (PBGs), analogous to the 
electronic band gaps in semiconductors, with in which the 
propagation of electromagnetic (EM) waves is completely 
forbidden. The existence of PBGs will lead to many in- 
teresting phenomena, e.g., modification of spontaneous 
emission [1-5] and photon localization [6-10]. Thus nu- 
merous applications of photonic crystals have been pro- 
posed in improving the performance of optoelectronic 
and microwave devices such as high-efficiency semicon- 
ductor lasers, right emitting diodes, wave guides, optical 
filters, high-Q resonators, antennas, frequency-selective 
surface, optical limiters and amplifiers [11-14]. These 
applications would be significantly enhanced if the band 
structure of the photonic crystal could be tuned. 

For the conventional PCs, the photonic band gaps re- 
main fixed once the PCs have been fabricated. If the 
band gaps of the photonic crystals could be tuned the 
applications would be significantly enhanced. A practical 
scheme for tuning the band gap was proposed by Busch 
and John [15]. It has been demonstrated theoretically 
and experimentally that PCs with liquid crystal infiltra- 
tion exhibit tunability on applying an external electric 
field [15]or changing the temperature [16-20]. If the con- 
stituent materials of PCs have magnetic permeabilities 
dependent on the external magnetic field, the photonic 
band gaps [PBGs] can be altered by changing the exter- 
nal magnetic field [21-28]. An electric or magnetic field 
changes the PBGs easily than the temperature. 

In Ref. [29], we have proposed special function pho- 
tonic crystals, which the medium refractive index is the 



function of space position, but the function value of re- 
fractive index is equal at two endpoints of every medium 
A and B, and obtain some new results different from 
the the conventional PCs. In this paper, we present a 
new general function photonic crystals (GFPCs), which 
refractive index is a arbitrary function of space position 
(needless refractive index same at two endpoint). Un- 
like conventional photonic crystals (PCs), which struc- 
ture grow from two materials, A and B, with different di- 
electric constants ea and sb- Firstly, we give the motion 
equation of light in one-dimensional GFPCs according to 
Fermat principle. Secondly, we calculate the transfer ma- 
trix for the one-dimensional GFPCs, which is different 
from the transfer matrix of the conventional PCs. Fi- 
nally, we give the dispersion relation, band gap structure 
and transmissivity. We choose the linearity refractive in- 
dex function for two medium A and B, and find the trans- 
missivity of GFPCs can be much larger or smaller than 1 
for the different slope linearity refractive index function, 
which is different from the transmissivity of conventional 
PCs (its transmissivity is in the range of and 1). Other- 
wise, we study the effect of different incident angles, the 
number of periods and optical thickness on the transmis- 
sivity, and obtain some new results. By the calculation, 
we find the conventional PCs is the special case of the 
GFPCs. 

2. The light motion equation in general func- 
tion photonic crystals 

For the general function photonic crystals, the medium 
refractive index is a periodic function of the space posi- 
tion, which can be written as n(z), n(x, z) and n(x, y, z) 
corresponding to one-dimensional, two-dimensional and 
three-dimensional function photonic crystals. In the fol- 
lowing, we shall deduce the light motion equations of the 
one-dimensional general function photonic crystals, i.e., 
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FIG. 1: The motion path of light in the medium of refractive 
index n(z). 



the refractive index function is n — n(z), meanwhile mo- 
tion path is on xz plane. The incident light wave strikes 
plane interface point A, the curves AB and BC are the 
path of incident and reflected light respectively, and they 
are shown in FIG. 1. 

The light motion equation can be obtained by Fermat 
principle, it is 



S [ n(z)ds = 0. 

J A 



(1) 



In the two-dimensional transmission space, the line ele- 
ment ds is 



ds = ^(dx) 2 + (dz) 2 = Vl + Pdx, (2) 
where z = , then Eq. (1) becomes 



S J n[z)y/l + (z) 2 dx = 0. 
The Eq. (3) change into 
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(3) 



'M»f*P)„ + »p fili . 0,(4) 



dz 



dz 



At the two end points A and B, their variation is zero, 
i.e., 5z(A) = 5z(B) = 0. For arbitrary variation Sz, the 
Eq. (4) becomes 



^l^TTT 2 -^z 2 (i + z 2 )-i 

dz dz 

, z\f\ + z 2 - z 2 'z(l + z 2 Y^ 



-n(z)- 



1 + z 2 

simplify Eq. (5), we have 

dn(z) zdz 



n(z) f + i 2 



0, (5) 



(6) 



The Eq. (6) is light motion equation in one-dimensional 
function photonic crystals. 

3. The transfer matrix of one-dimensional gen- 
eral function photonic crystals 
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FIG. 2: The light transmission and electric magnetic field 
distribution figure in FIG.l medium. 



In this section, we should calculate the transfer matrix 
of one-dimensional general function photonic crystals. In 
fact, there is the reflection and refraction of light at a 
plane surface of two media with different dielectric prop- 
erties. The dynamic properties of the electric field and 
magnetic field are contained in the boundary conditions: 
normal components of D and B are continuous; tangen- 
tial components of E and H are continuous. We consider 
the electric field perpendicular to the plane of incidence, 
and the coordinate system and symbols as shown in FIG. 
2. 

On the two sides of interface I, the tangential compo- 
nents of electric field E and magnetic field H are contin- 
uous, there are 



Eq — Ei — E t i + E' r2 . , 

H =H I = H n cos e{ - H' r2 cos B\. [ 1 



On the two sides of interface II, the tangential compo- 
nents of electric field E and magnetic field H are contin- 
uous, and give 

J En = E'j = Ei2 + E r2 

\H II = H' I = H t2 cos e\ - H r2 cos 0i, W 

the electric field E t \ is 

E n = E tw e^ XA+k ^ \ z=0 = E tw e 1 ^ sin 9 ' XA , (9) 

and the electric field E i2 is 

Ei2 — E tw e * * >\ z= b 

= E tw e i ^ n< - b ^ sineIiXB+cose * b \ (fO) 

Where xa and xb are x component coordinates corre- 
sponding to point A and point B. We should give the 
relation between Ei 2 and E t \. By integrating the two 
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FIG. 3: The structure (BA) N of the general function photonii 
crystals. 



sides of Eq. (6), we can obtain the coordinate compo 
nent of point B 



/•"(»> dn(z) _ f K * zdz 

Jn(0) n ( z ) Jko l + z 



2 " 



to get 



and 



dx = 



dz 



(11 



(12 



(13) 



where ko = cotO^ and k z = From Eq. (12), there is 
n(z) > n(0)sin9{. and the coordinate xb is 



Xb =x a + 



dz 



V / ( 1 + fc o)($}) 2 -l' 



(14) 



where 6 is the medium thickness of FIG. 1 and FIG. 2. 
By substituting Eqs. (9) and (14)into (10), and using the 
equality 



we have 



where 



n(0) sin0 t J = n(b) sin6>', 



E i2 = E n e iS \ 



5b = —n b (b)(cos9 
c 



Ib + smOl [ 
Jo 



dz 



(15) 



(16) 



= ), (17) 



ni sin 2 W 1 



and similarly 



E' = E r2 e^. 



(18) 



Substituting Eqs. (16) and (18) into (7) and (8), and 
using H = . I —nE, we obtain 
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FIG. 4: The two periods transmission figure of light in general 
function photonic crystals. 



where 



M E 



-in b (0). 



cos 5b 



cos d\ sin 5b 



i sin S b 



^j%n b (b) cos 61 . 

u(a)t =os el ™ a,_ |v>" 7 



n^, (£>) cos 1 . 



The Eq. (20) is the transfer matrix M in the medium of 
FIG. 1 and FIG. 2. By refraction law, we can obtain 



smtf, — 



n 



sin 0,° , cos 0/ = 
n(0) 



'1- 



n 2 (0) 



sin 2 ^, (21) 



where no is air refractive index, and n(0) = n(z)| z= o. 
Using Eqs. (15) and (21), we can calculate cos#f. 

4. The structure of one-dimensional general 
function photonic crystals 

In section 3, we obtain the M matrix of the half pe- 
riod. We know that the conventional photonic crystals 
is constituted by two different refractive index medium, 
and the refractive indexes are not continuous on the in- 
terface of the two mediums. We could devise the one- 
dimensional general function photonic crystals structure 
as follows: in the first half period, the refractive index 
distributing function of medium B is rib(z). and in the 
second half period, the refractive index distributing func- 
tion of medium A is n a (z), corresponding thicknesses are 
b and a, respectively. Their refractive indexes satisfy con- 
dition rib(b) 7^ n a (0), their structure are shown in FIG. 
3, and FIG. 4. The Eq. (20) is the half period transfer 
matrix of medium B. Obviously, the half period transfer 
matrix of medium A is 



cos 5 a 



i sin S a 



M A = 



-m o (0) 



^ cos Ol 1 sin S a 



n a (a) cos Of* 



,(0) 
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where 



Sa 



— n a {a)[cos6 I i I ■ a 

+sine ii r dz = ] 

Jo I i 



0?, 



and 



sin^ 7 = ^-sin^ , 
n (o) 



cos0/ J = Wl- 



„:„2 ^0 



■ sin 



n2(a) 

In one period, the transfer matrix M is 
M = M B ■ M A 



cosS b -' sin<5 " ■ 

/^n»(6)co.«f 



) COS 0] 

— 2 sin (5 a 



008 "7S^ 

V w v 

-tn„(O),/^cos0"sin<y o "° (0) . cos9 ' 



cos , 
> t " cos (5, 



t) cos } 



(77 



(23) 



(24) 



(25) 



(26) 




The form of the GFPCs transfer matrix M is more com- 
plex than the conventional PCs. The angle 0{, 9f , Oj. 1 
and Of 1 are shown in Fig. 4. The characteristic equation 
of GFPCs is 



= M 1 M 2 • • • M 



AT 



En+i 
Hn+i 



= M b M a M b M a ■ ■ ■ M b M a ^ JV+i 
En+i \ _ f A B \ / 



En+i 



M 



JV+l 



H 



N+l 



(28) 



Where TV is the number of period. 

5. The dispersion relation, band gap structure 
and transmissivity 



With the transfer matrix M (Eq. (27)), we can study 
the dispersion relation and band gap structure of the 
function photonic crystals. 
Since 



En ] = M ( En+1 
Hn J \ Hn+i 



By Bloch law, there is 
En 



\/„.U„( ). (29) 



-ikd 



En+i 
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(30) 
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FIG. 5: The line refractive index function in a period. The 
straight line (a-e) are corresponding to different line functions 
for medium A. 



where d = b + a. With Eqs. (29) and (30), there is 



^ l-M 6 M a l ^ + ]-e I ),(31) 



The non-zero solution condition of Eq. (31) is 

det(M b M a - e- lkd ) = 0, (32) 

i.e., 



(cos 5 b cos 5 a sin 5 b sin 5 a — e 1 ) (cos 5 b cos 5 a sin 5 b sin 5 a — e 1 ) 

Vb Va 

% % 

+ ( cos<5fc sin<5 a smS b cos 5 a )(— ir) b sin5 b cos<5 a — ir\ a cos5 b sin5 a ) = 0. (33) 

Va Vb 



I 

Simplifying Eq.(33), we obtain the dispersion relation From Eq. (34), we can study the photonic dispersion re- 

cos fed = cos 5 b cos 5 a — -(— + — ) sin sin 5 a . (34) 

2 Va Vb 




FIG. 6: Comparing the transmissivity of the GFPCs with 
different line refractive index function: Fig. 6 (a-e) corre- 
sponding to Fig. 5 (a-e). 



FIG. 7: Comparing the transmissivity of the GFPCs with 
different line refractive index function for medium B. 



lation and band gap structure, and can obtain the trans- 
mission coefficient t from Eq. (28) 

t = EtN + 1 = ^ (-35) 

En Ai] + BrjoriN+i + C + Drj N+1 ' 

and transmissivity T is 

T = t-t*. (36) 

6. Numerical result 

In this section, we report our numerical results of trans- 
missivity. We consider refractive indexes of the linearity 
functions in a period, it is 

n b (z)=n b (0) + nb{b) ~ nb{0) z, 0<z<6, (37) 
o 

n a (z) = n a (0) + na{a) - na{0) z, 0<z<a, (38) 

Eqs. (37) and (38) are the refractive indexes distribution 
function of two half period mediums B and A , which 
are shown in FIG. 5. When the endpoint values n b (Q), 
n-b(b), n a (0) and n a (a) are all given, the line refractive 
index functions n b (z) and n a (z) are ascertained. 

The main parameters are: the half period thickness 
b and a, the starting point refractive indexes n b (0) and 
n a (0), and end point refractive indexes n b (b) and n a (a), 
the optical thickness of the two mediums are equal, i.e., 
n b (0)b = n a (0)a. 

In Fig. 6, we take n b (0) — 1.38, n b (b) — 1.9, i.e., the 
refractive index distribution of medium B is conformable, 
n a (0) = 2.35, the incident angle 8® = 0, the center fre- 
quency loq = 1.215 x 10 15 Hz, A = the thickness 
b = 280nm and a — 165nm and the period number 



N = 16. From Fig. 6(a-e), the end point value of n a (a) 
are taken as: 4.25, 2.6, 2.35, 2.1 and 1.5 respectively. 
And the corresponding refractive indexes distribution are 
straight lines (a), (b), (c), (d) and (e) in Fig. 5. When 
the refractive index distribution of medium B maintain 
unchanging, with the end point value n a (a) of medium A 
decrease, we can obtain the following results for the one- 
dimensional GFPCs: (1) The maximum of transmissiv- 
ity decrease. From Fig. 6 (a) to (e), their transmissivity 
maximums are 2500000, 900, 180, 30 and 0.15, respec- 
tively. (2) The width of band gaps become narrow and 
the number of band gaps decrease. (3)The transmissivity 
of one-dimensional GFPCs can be much larger or smaller 
than 1 for different slope linearity refractive index func- 
tion, which is different from the transmissivity of con- 
ventional PCs (its transmissivity T < 1). When we take 
n b = n b (b) and n a = n a (a), and calculate the trans- 
missivity by the GFPCs equations, i.e., Eqs. (27) and 
(35), the calculate results of transmissivity is the same 
as the conventional PCs. So, the conventional PCs is the 
special case of the GFPCs. 

In Fig. 7, we take n a (0) = 2.35, n a (a) = 2.4, i.e., the 
refractive index distribution of medium A is confirmablc, 
rifc(O) = 1.38, and the end point value n b (b) of medium B 
are taken as: 1.9 and 1.5. By the two refractive indexes 
distributions of medium B, we obtain the transmissiv- 
ity are corresponding to Fig. 7(a) and (b). It can be 
found that when the refractive index of medium A keep 
unchanging, with the end point value n b (b) of medium 
B decrease, the transmissivity maximum and the num- 
ber of band gaps decrease, and the width of band gaps 
increase. 

In Fig. 8, Fig. 9 and Fig. 10, the main parameters are: 
n h (0) = 1.38, n h (b) = 1.9, n a (0) = 2.35, n a (a) = 2.6. In 
the inset of Fig. 8, we discuss the influence of period 
number N on the transmissivity. The results shown in 
Fig. 8 (a-c) are calculated by the period number N = 12, 
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N = 16 and N — 20, respectively. We can see that when 
the period number of the FGPCs is small, for example, 
N = 12 the maximum of transmission intensity achieve 
160 [Fig. 8(a)]. When the period number increase up to 
20 [Fig. 8(c)], the transmission intensity is nearly 5000, 
it is much higher than N = 12, i.e., as the period number 
increase, the transmission intensity increase. 

Figure 9 shows that the effect of different incident an- 
gles on transmission intensity. In Fig. 9 (a-d), the inci- 
dent angles are: 0? = 0, 0? = § , 0? = f and 0? = § , 
respectively. As we can see in Fig. 9 (a-d), when the in- 
cident angle increase, the transmission intensity increase, 
and the number of band gaps increase. 

Figure 10 shows the variation of transmissivity due to 
the change of the optical thickness. The optical thickness 
of A and B medium are equal, i.e., n&(0)6 = n a (0)a. 
The results shown in Fig. 10 (a-c) are calculated by the 
optical thickness as: ^ and 4f- respectively. It is 
found when the optical thickness increase, the number of 
band gaps increase, while the maximum of transmission 
intensity keep unchanged. 

In the following, we consider two groups line refractive 
indexes functions in a period, which are all decreasing 
linearly, they are shown in Fig. 11. The first group main 
parameters are: n&(0) = 1.9, nt,(b) = 1.8, n a (0) = 2.6 
and n a (a) = 2.5, it is the solid line (a). The second 
group main parameters are: nb(0) = 1.5, nt(b) = 1.4, 
n a (0) = 2.4 and n a (a) — 2.3, it is the dash-dot line (b). 
The optical thickness are same: n&(0)6 = n a (0)a = 

Figure 12 is the transmission spectra. The Fig. 12 (a) 
correspond to the first group line refractive indexes func- 
tions, the Fig. 12 (b) correspond to the second group line 
refractive indexes function. Here, we notice that when 
the refractive indexes function of the two mediums A and 
B are all decreasing, the transmission intensity decrease 
and the width of band gaps become wider. 

In Fig. 13-15, the main parameters are: n&(0) = 1.9, 
n b (b) = 1.8, n o (0) = 2.6 and n a (a) = 2.5. 

In Fig. 13, we study the effect of period number TV 
on the transmission intensity. The incident angle Q\ = 0, 
the optical thickness nb(0)b — n a (0)a = They are 
shown in Fig. 13 (a-c), and the period number are N = 
12, N — 16 and N — 20, respectively. It can be found 
that as period number N increases from 12 to 20, the 
transmission intensity decrease. 

In addition, we study the transmission intensity with 
different incident angles as shown in Fig. 14. We 
take the period number N — 16, and optical thickness 
nb(0)b = n a (0)a = we can find when the incident an- 
gle 0° increase, the transmission intensity decrease. For 
example, when the incident angle 0° = 0, the maximum 
of transmission intensity nearly 0.08 (Fig. 14(a)), and 
0? = -|, the maximum of transmission intensity nearly 
0.08 (Fig. 14(d)). 

In what follows, we will take into account the trans- 
mission at different optical thickness, which is shown in 
Fig. 15. We also choose 0? = 0, N = 16. When the 
optical thickness are taken ^ and 4p, the transmis- 
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FIG. 8: Comparing the transmissivity of the function PCs 
with different period number iV: (a) N=12, (b) N=16 and (c) 
N=20. 
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FIG. 9: Comparing the transmissivity of the GFPCs with 
different incidence angle: (a)0? = 0, (b) 0° = f , (c)0? = f 
and (d) e°i = f . 



sion spectra are corresponding to Fig. 15 (a-c). we can 
find when the optical thickness increases, the number of 
band gaps increase, the width of band gaps become nar- 
row, while the maximum of transmission intensity keep 
unchanged. 

7. Conclusion 

In summary, We have theoretically investigated a new 
kind of general function photonic crystals (GFPCs), 
which refractive index is a function of space position. 
Based on Fermat principle, we achieve the motion 
equations of light in one-dimensional general function 
photonic crystals, and calculate its transfer matrix. We 
choose the linearity refractive index function for two 
mediums A and B, and find the transmissivity of one- 
dimensional GFPCs can be much larger or smaller than 
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FIG. 10: Comparing the transmissivity of the GFPCs with 
different optical thickness: (a)^, (b) ?f an d ( c )^. 



FIG. 12: Comparing the transmissivity of the GFPCs with 
different line refractive index function: (a) the first group 
function and (b) the second group function. 
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FIG. 11: The two groups line refractive index function in a 
period: (a) the first group function and (b) the second group 
function. 



1 for different slope linearity refractive index function. 
Otherwise, we study the effect of different incident an- 
gles, the number of periods and optical thickness on the 
transmissivity, and obtain some new results as follows: 

(1) when nt,(b) > rib(0) and n a [a) > n a (0), the peak 
value of transmission intensity is much larger than 1. 

(2) when n^b) < rib(0), n a (a) < n a (0), the transmission 
intensity is much smaller than 1. With the high and low 
transmissivity of GFPCs, we can make amplification and 
decay device of light. (3) For arbitrary medium A or B, 
when its slope of line refractive index function decrease 
the transmissivity peak value decrease, the number of 
band gaps decrease and the band gaps become narrow. 
(4) When the transmissivity T > 1, as the optical 
thickness increase, the number of band gaps increase, 
and its width become narrow, while the peak value of 
transmissivity is unchanged. (5) When the transmissiv- 



FIG. 13: Comparing the transmissivity of the function PCs 
with different periodicity (a) N=12 (b) N=16 (c) N=20 
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FIG. 14: Comparing the transmissivity of GFPCs with dif- 
ferent incidence angles: (a)6»° = 0, (b) 8° = f , (c)0? = f and 
(d) 0° = §. 



8 




0.2 



0.1 



0.0 



llll UfatiJUli milto j|g|Hly ||gjJL iiglll 

^nnnrnrnrnrn 



(c) 4 



'08 



FIG. 15: Comparing the transmissivity of the GFPCs with 
different optical thickness: (a)^, (b) ^ and (c)^. 



ity T < 1, as the optical thickness increase, the number 
of band gaps increase, and its width become narrow, 
while the peak value of transmissivity is unchanged. (6) 
When the transmissivity T > 1, as the incident angles 
increase, the peak value of transmissivity also increase, 
and the position of band gaps red shift. (7) When the 
transmissivity T < 1, as the incident angles increase, the 
peak value of transmissivity decrease, and the position 
of band gaps red shift. (8) When the transmissivity 
T > 1, as the number of period increase, the peak value 
of transmissivity increase, while the position and width of 
band gaps are unchanged. (9) When the transmissivity 
T < 1, as the number of period increase, the peak value 
of transmissivity decrease, while the position and width 
of band gaps are unchanged. The new results of GFPCs 
couldn't be found in conventional PCs, and we think the 
GFPCs should be widely applied in the future. 
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